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The Standard Approach (SA) for description of the structure function gi combines the DGLAP 
evolution equations and Standard Fits for the initial parton densities. The DGLAP equations 
describe the region of large and large x, so there are not theoretical grounds to exploit them at 
small X. In practice, extrapolation of DGLAP into the region of large and small x is done with 
complementing DGLAP with special, singular (~ a;~") phenomenological fits for the initial parton 
densities. The factors x''^ axe wrongly believed to be of the non-perturbative origin. Actually, they 
mimic the resummation of logs of x and should be expelled from the fits when the resummation is 
accounted for. Contrary to SA, the resummaton of logarithms of a; is a straightforward and natural 
way to describe gi in the small-x region. This approach can be used at both large and small 
r-j , where DGLAP cannot be used by definition. Confronting this approach and SA demonstrates that 

^ the singular initial parton densities and the power Q^-corrections (or at least a sizable part of them) 

I ' are rather not real physical phenomena but the artefacts caused by extrapolating DGLAP into the 

small-x region. 

r^, ' PACS numbers: 12.38.Cy 

> ■ I. INTRODUCTION 

The Standard Approach (SA) for description of the structure function gi involves the DGLAP evolution equations[l| 
and Standard Fits 2] for the initial parton densities 5q and 5g. The fits are defined from phenomenological consider- 
ations at a; ~ 1 and = fi^ ^ IGeV^. The DGLAP equations are one-dimensional, they describe the -evolution 
only, converting 6q and Sg into the evolved distributions Aq and Ag. They represent gi at the region A: 



o . 

^ ■ A: Q2 ^ ^ < ^_ 

The X -evolution is supposed to come from convoluting Aq and Ag vi^ith the coefficient functions Cdglap- However, 
in the leading order C^'qj^^p = 1; the NLO corrections account for one- or two- loop contributions and neglect higher 
^\ loops. It is the correct approximation in the region A but becomes false in the region B: 

H ; 

. F^. . B: Q2^^2^ x<l (2) 

where contributions ~ ln'^(l/a;) are large and should be accounted for to all orders in as- Cdglap do no include 
the total resummation of leading logarithms of x (LL), so there are not theoretical grounds to exploit DGLAP at 
small X. However regardless of that, SA extrapolates DGLAP into the region B, invoking special fits for 6q and 6g. 
A general structure of such fits (see Refs. 2]) is as follows: 

Sq = Nx-^ifix) (3) 

where iV is a normalization constant; a > 0, so x^"' is singular when x — > and (p{x) is regular in a; at a; —> 0. In 
Ref. Q we showed that the role of the factor x^"^ in Eq. ^ is to mimic accounting for the total resummation of LL 
performed in Refs [4|, jsj. Similarly to LL, the factor x~"' provides the steep rise to gi at small x and sets the Regge 
asymptotics for gi at x —> 0, with the exponent a being the intercept. The presence of this factor is very important 
for extrapolating DGLAP into the region B: When the factor a;~° is dropped from Eq. ([3]), DGLAP stops to work at 
X < 0.05 (see Ref. 3 fo^' detail). Accounting for the LL resummation is beyond the DGLAP framework because LL 
come the phase space violating the base of DGLAP: the DGLAP -ordering 

^I^<kl^<kl^< ... < q2 (4) 
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for the ladder partons. LL can be accounted only when the ordering Eq. (|4]) is lifted and all ki j_ obey 

<kf^<{p + qf w (1 - x)2pq « 2pq 
at small x. Replacing Eq. ([4]) by Eq. ([5]) leads inevitably to the change of the DGLAP parametrization 



(5) 



(6) 



by the alternative parametrization of as given by Eq. (jl3p . This parametrization was obtained in Ref. @ and was 
used in Refs. d, Q in order to find explicit expressions accounting for the LL resummation for gi in the region B. 
Obviously, those expressions invoke the fits for the initial parton densities without the singular factors a;~". Let us 
note that replacement of Eq. @ by Eq. ([5]) brings a more involved /i -dependence to gi. Indeed, Eq. @ makes 
contributions of gluon ladder rungs be infrared (IR) stable, with ^ acting as a IR cut-off for the lowest rung and ki ± 
playing the role of the IR cut-off for the i + 1-rung. In contrast, Eq. ([5]) implies that ji acts as the IR cut-off for every 
rung. 

Besides the regions A and B, it i necessary to know gi in the region C: 



(7) 



because this region is studied experimentally by the COMPASS collaboration. Obviously, DGLAP cannot be exploited 
here. Alternatively, in Refs. [3, Q we obtained expressions for gi in the region C. In Ref. 7] we showed that gi 
practically does not depend on x at small x, even at x ^ 1. Instead, it depends on the total invariant energy 2pq. 
Experimental investigation of this dependence is extremely interesting because according to our results gi, being 
positive at small 2pq, can turn negative at greater values of this variable. The position of the turning point is sensitive 
to the ratio between the initial quark and gluon densities, so its experimental detection would enable to estimate 
this ratio. In Ref. Q we analyzed the power contributions ~ to gi usually attributed to higher twists. We 

proved that a great amount of those corrections have a simple perturbative origin and resummed them. Therefore, 
the genuine impact of higher twists can can be estimated only after accounting for the perturbative -corrections. 



II. DESCRIPTION OF gi IN THE REGION B 

The total resummation of the double-logarithms (DL) and single- logarithms of x in the region B was done in 
Refs. [1, In particular, the non-singlet component, g^^ of gi is 



g^^'ix, Q') = (e^/2) / (l/x)-C7^5(^)59(t^) exp (i/jvs(t^) HQVf^^)) , 



(8) 



with new coefficient functions Cats, 



and anomalous dimensions H 



NS, 



(9) 



where 



Hns = (1/2) w - - B{io) 



B{lj) = (47rCF(l +tj/2)A(cj) -|-D(w))/(27r2) 



(10) 



(11) 



D{u}) and A{lo) in Eq. (fTT|) are expressed in terms of p = ln(l/.T), rj = ln(/i^/Agp^), b = (33 — 2n/)/127r and the 
color factors Cp = 4/3, = 3: 



p + ri i 
{p + 77)^ + TT^ 77 



(12) 



Ai.) = - 



if + TT^ [p + ry)2 + 



(13) 
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Hs and Cms account for DL and SL contributions to all orders in as- Eq. ((T3)) and ([12]) depend on the IR cut-off /i 
through variable 77. It is shown in Refs. [1, Q that there exists an Optimal scale for fixing /x: /x « 1 Gev for g^^ and 
/i ~ 5 GeV for gl . The arguments in favor of existence of the Optimal scale were given in Ref. Q . Eq. ([8]) predicts 
that gi exhibits the power behavior in x and when a; — > 0: 

gr - {QVxY^'/\ 5f - {QVx^^/' (14) 

where the non-singlet and singlet intercepts are A^s = 0.42, Ag = 0.86 respectively. However the asymptotic 
expressions should be used with great care: According to Ref. [3], Eq. should not be used at a; > 10~^. So, 
Eq. ([8]) should be used instead of Eq. ([14]) at available small x. Expressions accounting the total resummation of LL 
for the singlet gi in the region B were obtained in Ref. [H]- They are more complicated than Eq. ([8]) because involve 
two coefficient functions and four anomalous dimensions. 



III. DESCRIPTION OF gi IN THE REGION C 

Region C is defined in Eq. ([7]). It includes small Q^, so there are not large contributions ln'"'((5^//z^) in this region. 
In other words, the DGLAP ordering of Eq. ^ does not make sense in the region C , which makes impossible 
exploiting DGLAP here. In contrast, Eq. ^ is not sensitive to the value of and therefore the total resummation 
of LL does make sense in the region C. In Ref. Q we suggested that the shift 

Q^^Q^+fi^ (15) 

would allow for extrapolating our previous results (obtained in Refs. [1, [E] for <7i in the region B) into the region C. 
Then in Ref. [1] we proved this suggestion. Therefore, applying Eq. (|15p to g^^ leads to the following expression for 
gi^ valid in the regions B and C: 

g^^{x + z,Q') = {el/2) ^ ^ (_l_)"cArs(c.)<5q(c.) exp (ff^s(c.) In ((Q^ + /i^)^')) , (16) 

where z = /2pq. Obviously, Eq. ([16]) reproduces Eq. ([8]) in the region B. Expression for gf looks similarly but more 
complicated, see Refs. [3, [81] for detail. Let us notice that the idea of considering DIS in the small-Q^ region through 
the shift Eq. (|15p is not new. It was introduced by Nachtmann in Ref. 10] and used after that by many authors (see 
e-g- [HI) I being based on different phenomenological considerations. On the contrary, our approach is based on the 
analysis of the Feynman graphs contributing to gi . 



IV. PREDICTION FOR THE COMPASS EXPERIMENTS 



The COMPASS collaboration now measures the sing let gf sX X ^ 10^3 and < 1 GeV^ i.e. in the kinematic 
region beyond the reach of DGLAP. However, our formulae for g'^^ and gf obtained in Refs. [3, Q cover this region. 
Although expressions for singlet and non-singlet gi are different, with formulae for the singlet being much more 
complicated, we can explain the essence of our approach, using Eq. (|16p as an illustration. According to results of 0], 
/X « 5 GeV for gf, so in the COMPASS experiment <^ /i^. It means, \n~{Q^ + fi"^) can be expanded into series in 
Q'^/lJ?, with the first term independent of Q^: 

gfix + z, Q^ ^J?) = gf{z, ^?) + Y.{Q^ / Ek{z) (17) 
where Ek{z) account for the total resummation of LL of z and 

/too 7 
_ ^^i^/'T[cU^)^ii^) + cu^)H^)]^ (18) 

so that 5q{u!) and 5g{uj) are the initial quark and gluon densities respectively and C|'^ are the singlet coefficient 
functions. Explicit expressions for Cg^ are given in Refs. The standard fits for Sq and 6g contain singular 

factors ^ x~°' which mimic the total resummation of leading logarithms of x. Such a resummation leads to the 
expressions for the coefficient functions different from the DGLAP ones. After that the singular factors in the fits can 
be dropped and the initial parton densities can be approximated by constants: 

Sq « iV, Sg^Ng, (19) 
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FIG. 1: Gi evolution with decreasing z = /2[pq) for different values of ratio r — Sg/5q: curve 1 - for r = 0, curve 2 - for 
r — —5 , curve 3 -for r = —8 and curve 4 -for r = — 15. 

SO, one can write 

51 (Q' < m') ~i<el> /2)Nfii{z) (20) 

with 

— (l/zr[C| + (7V,/7V,C|)]. (21) 

Obviously, Gi depends on the ratio Ng/Nq. The results for different values of the ratio r = Ng/Nq, Gi are plotted 
in Fig.[T] When the gluon density is neglected, i.e. Ng = (curve 1), Gi being positive at x ^ 1, is getting negative 
very soon, at 2; < 0.5 and falls fast with decreasing z. When Ng/Nq — —5 (curve 2), Gi remains positive and 
not large until z ^ 10^^, turns negative at z ~ 0.03 and falls afterwards rapidly with decreasing z . This turning 
point where Gi changes its sign is very sensitive to the magnitude of the ratio r . For instance, at Ng/Nq = —8 
(curve 3), Gi passes through zero at 2: ^ 10"'^. When Ng/Nq < —10, Gi is positive at any experimentally reachable 
z (curve 4) . Therefore, the experimental measurement of the turning point would allow to draw conclusions on the 
interplay between the initial quark and gluon densities. 

V. REMARK ON THE HIGHER TWISTS CONTRIBUTIONS 

In the region B one can expand terms ^ {Q^ + /i^)*^ in Eq. (fT6|) into series in (/i^/Q^)" and represent g^^{x + 
z,Q'^,y/^) as follows: 

g^ix + z,Q\ = g^H^, QVa^') + Y^^^f^'/Q'^n (22) 

fc=i 

where g^ ^ {x , Q'^ / /s^) is given by Eq. ([5]); for explicit expressions for the factors Tk see Rcf. The power terms in 
the rhs of Eq. look like the power ^ 1/(Q^)'^ -corrections and therefore the Ihs of Eq. ([22)1 can be interpreted 
as the total resummation of such corrections. These corrections are of the perturbative origin and have nothing in 
common with higher twists contributions (= HTW). The latter appear in the conventional analysis of experimental 
date on the Polarized DIS as a discrepancy between the data and the theoretical predictions, with g^ ^ {x , / ^'^) 
being given by the Standard Approach: 

gNS exp ^ ^NSSA ^ ^rj.^ (23) 

Confronting Eq. to Eq. leads to an obvious conclusion: For estimating genuine higher twists contributions 
to gi^ , one should account, in the first place, for the perturbative power corrections predicted by Eq. (j^ : otherwise 
the estimates cannot be reliable. It is worth mentioning that we can easily explain the empirical observation made 
in the conventional analysis of experimental data: The power corrections exist for > 1 GeV^ and disappear when 
^ 1 GeV^. Indeed, in Eq. (gH) ^ = 1 GeV , so the expansion in the rhs of Eq. make sense for > 1 GeV^ 
only; at lesser it should be replaced by the expansion of Eq. ^TE\i in {Q'^/fi^y\ 
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VI. CONCLUSION 

Resummation of the leading logarithms of x is the straightforward and most natural way to describe gi at small 
X. Contrary to DGLAP, our approach is not sensitive to the value of and allows one to describe gi at small x 
and arbitrary in terms of the same expressions at large and small Q^. We have used it for studying the gi singlet 
at small because this kinematic is presently investigated by the COMPASS collaboration. It turns out that gi in 
this region depends on z = fi^ /2pq only and practically does not depend on x, even at x <C 1. Numerical calculations 
show that the sign of gi is positive at z close to 1 and can remain positive or become negative at smaller z, depending 
on the ratio between Sg and 5q. It is plotted in Fig. 1 for different values of Sg/Sq. Fig. 1 demonstrates that the 
position of the sign change point is sensitive to the ratio Sg/6q, so the experimental measurement of this point would 
enable to estimate the impact of 6g. 

The alternative to the resummation is extrapolating DGLAP from its natural region of applicability (large x and 
large Q^) into the region of small x and large Q^. As the DGLAP equations cannot account for the LL resummation, 
SA mimics the resummation through the special choice of the fits for the initial parton densities: the singular factors 
in the fits cause the steep rise of gi at small x and provide the Regge asymptotics for gi (however with the incorrect 
phenomenological intercepts) when x — > 0. They should be dropped when the total resummation of LL of x is taken 
into account. The remaining, regular x-terms of the DGLAP fits (the function in Eq. ([3])) can obviously be replaced 
by much simpler expressions, so the number of phenomenological parameters in the fits can be reduced from 5 to 2 or 
even 1. To conclude, let us notice that extrapolating DGLAP into the small-x region, though provides a satisfactory 
agreement with experimental data, leads to various wrong statements. We enlisted the most of them in a recent 
Ref. "9*1. Below we mention two more such wrong statements: 

Statement 1: The -power corrections stem from higher twists gi and can he measured as the discrepancy between 
the DGLAP predictions and the data. This statement is wrong as shown in the previous Sect. 

Statement 2: The impact of the LL resummation on the small-x behavior of gi is small. This statement appears 
when the resummation has been included into the DGLAP expressions where the fits contain singular factors. Such 
inclusion is inconsistent and means actually a double counting of the LL contributions: once through the fits and 
secondly in the explicit way. It also affects the small-x asymptotics of gi, leading to the incorrect values of the 
intercepts of gi (see Ref. Q for more detail) . 
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